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ABSTRACT 
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Oh 

I Context. The Cosmic Far-Infrared Background (CIB) at wavelengths around 160 jim corresponds to the peak intensity of the whole 
Extragalactic Background Light, which is being measured with increasing accuracy. However, the build up of the CIB emission as a func- 
tion of redshift, is still not well known. 

Aims. Our goal is to measure the CIB history at 70 yum and 160 yum at different redshifts, and provide constraints for infrared galaxy evolution 
models. 

Methods. We use complete deep Spitzer 24 fim catalogs down to about 80 /jJy, with spectroscopic and photometric redshifts identifications, 
from the GOODS and COSMOS deep infrared surveys covering 2 square degrees total. After cleaning the Spitzer/MIPS 70 fjm and 160 /im 
maps from detected sources, we stacked the far-IR images at the positions of the 24 /jm sources in different redshift bins. We measured the 
contribution of each stacked source to the total 70 and 160 /jm light, and compare with model predictions and recent far-IR measurements made 
■ ' with Herschel/PACS on smaller fields. 
f ^ , Results. We have detected components of the 70 and 160 yum backgrounds in different redshift bins up to z ~ 2. The contribution to the CIB is 
maximum at 0.3 < z < 0.9 at 160yum (and z < 0.5 at 70 fim). A total of 81% (74%) of the 70 (160) pm background was emitted at z < 1. We 
estimate that the AGN relative contribution to the far-IR CIB is less than about 10% at z < 1.5. We provide a comprehensive view of the CIB 
buildup at 24, 70, 100, 160 fim. 

Conclusions. IR galaxy models predicting a major contribution to the CIB at z < 1 are in agreement with our measurements. The consistency of 
our results with those obtained through the direct study of Herschel far-IR data at 160 fim confirms that the stacking analysis method is a valid 
approach to estimate the components of the far-IR background using prior information on resolved mid-IR sources. Our results are available 
online http://www.ias.u-psud.fr/irgalaxies/ . [(...) abstract abridged] 
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C$ ' Key words. Cosmology: observations, Diffuse Radiation - Galaxies: Evolution, Starburst, Active Galactic Nuclei, Infrared, BL Lacertae 
objects 



1. Introduction 

The Extragalactic Background Light (EBL) is the relic emis- 
sion of galaxy formation and evolution, i.e. due to star forma- 
tion and accretion processes (with this definition, the Cosmic 
Microwave Background due to recombination at redshift z ~ 
1100 is not part of the EBL). The EBL spectrum peaks 
in the Far-Infrared (FIR), where it is commonly known as 
the Cosmic Infrared Background (CIB) ( [Puget et al., 1 996 
IHauser et al.,T998l |Hauser & Dwek, 200 1[ |Kashlinsky, 2 005 
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|Dole et al., 2006| l. The EBL and the CIB encode the emission 
processes of structure formation, and can thus be used to con- 
strain the photon budget of the cooling processes leading the 
baryons to fall within the dark matter halos and form galaxies. 
The measurements of the EBL level and structure bring thus 
one of the many useful constraints for the models. 

The CIB Spectral Energy Distribution is known with 
increasing accuracy (for instance in the FIR and submil- 
limetre regime: Puget et al., 1996 |Aharonian et al., 2006 



|Dole et al., 2006| [Bethermi n et al., 2010] >, but little is known 
about its history, i.e. its buildup as a function of redshift. 
This missing information should help constrain galaxy 
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Fig. 1. Number counts at 24 /an in the GOODS HDFN (blue 
diamond), GOODS CDFS (gray triangle) and COSMOS (pink 
square) fields. The errors bars used only include Poisson statis- 
tics, and not cosmic variance. 
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Fig. 2. Cumulative stacked brightness at 70 /mi and 160 /im 
(in nW.nr 2 .sr~') on the CLEANed maps, as a function of the 
24 /an flux of our sample, regardless of the redshift. the 3 
fields are represented at 160 and 70 /an: COSMOS; GOODS-N 
(HDFN) and GOODS-S (CDFS). 



evolution models, and also better understand the physics 
of blazars, whose high-energy photons interact with the 
CIB along the line of sight (e.g. |Aharonian et al., 2 007 
Albert & Magic Collabo ration, 2008| |Raue et al., 2009] 

Kneiske & Dole, 2009) . 



The history of the CIB buildup can be derived by in- 
tegrating the luminosity functions of galaxies as a function 
of redshift (neglecting other sources of diffuse emission and 
thus assuming that the CIB is due to galaxies). This is a 
very difficult task in practice, since high-redshift luminos- 
ity functions are not yet measured at wavelengths close to 
the peak of the CIB (near 160 /mi) but instead in the mid- 
infrared range (e.g. |Le Floc'h et al., 20 05 Capu ti et al., 2007) , 
or only in the local universe ( Soifer & Neugebauer, 1991| 
Take uchi et al., 2006) . This situation is about to change 
with the latest Spitzer surveys and the ongoing deeper 
Herschel surveys ( |Magnelli et al., 2009| |Clements et al., 2010| 
Dye et al., 2010) . 



Two recent breakthroughs have been made by using 
COSMOS and GOODS surveys. Firstly, using about 30000 
Spitzer 24 /mi selected sources with accurate photometric red- 
shifts dllbert et al„ 2009) , [LeFToc'h et al. (2009)1 were able to 
measure the 24 //m background buildup with redshift (e.g. 
their figures 7 to 9). They furthermore show that the redshift 
information is crucial when confronting data to the models, 
since it helps breaking degeneracies. Secondly, using the red- 
shift identification of Herschel/PACS 100 and 160 /mi sources, 
|Berta et al. (2010)] were able to measure the CIB build up in 
four redshift bins, in the 140 arcmin 2 GOODS-N field (an area 
about 40 times smaller than used in this analysis). 



In this paper, we measure the 70 /mi and 160 /mi CIB his- 
tory since z = 2, with the use of a stacking analysis of galaxies 
detected at 24 //m (a good proxy for the 160 //m CIB popu- 
lation, e.g. |Dole et al. (2006)| |Bethermin et al. (2010)) in the 
Spitzer data of the GOODS and COSMOS fields. This ap- 
proach complements on a large area what is being done with 
Herschel in |Berta et al. (2010)| at 100 and 160 /an. 

2. Data and Sample 

2.1. GOODS data 

The data were acquired by the MIPS imaging photometer 
at 24 /mi, 70 /mi and 160 /mi ( [Rieke et al., 2 004) onboard 
the Spitzer infrared space telescope ( Werner et al., 2004), 
and come from the GOODS tea m (|Chary et al., 2004 
and Guaranteed Time observations (Papovich et al., 2004 
|Dole et al., 2"004) of the Chandra Deep Field South 
(CDFS) and the Hubble Deep Field North (HDFN). 



Papovich et al. (2004) extracted a catalog at 24 /an, with 



80% completeness at 80 /iJy. We use a sample of 1349 galaxies 
with 24 fim flux densities S 24 > 80 /Jy, located in the two 
GOODS fields, north and south, for a total area of 291 Sq. 



Arcmin ( |Caputi et al., 2006 Caputi et al., 2007). The galaxies 
have been completely identified, and redshifts determined 
for all of them, with more than 45% of spectroscopic red- 
shifts. Active Galactic Nuclei (AGN) are separated from the 
star-forming systems using X-ray data and near-infrared (3.6 
to 8 /mi) colors: we have 136 AGNs for 1213 star-forming 



systems ( [Caputi et al., 2007 1. For our purpose of measuring 
the contribution of mid-infrared galaxies to the far-infrared 
background by redshift slice, we cut the 24 /mi sample in four 
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Fig. 3. Images of all the stacked galaxies on the 160 fim (top) and 70 fim (bottom) GOODS CLEANed maps by redshift bin (right 
to left): 0<z<0.65; 0.65<z<1.3; 1.3<z<2 and z>2. The S/N ratio is indicated in each image. Notice the detection in the two first 
redshift bins at both wavelengths. Images are 305 x 305 sq. arcsec. wide at 70 fim (with 9.85 arcsec pixel plate), and 496 x 496 
sq. arcsec. at 160 /mi (with 16 arcsec pixel plate). The PSF FWHM being 18 arcsec. (resp. 40) at 70 fim (resp. 160 fim), the PSF 
shown on these figures has about the same extend of 1.8 to 2.5 pixels at both wavelengths. 
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Fig. 4. Images of all the stacked galaxies on the 160 fim COSMOS CLEANed maps with 14 redshift bins (left to right and top 
to bottom). The S/N ratio is indicated in each image. Notice the detection in all redshift bins at both wavelengths. Images are 
488 x 488 sq. arcsec. at 160 fim (with 8 arcsec pixel plate). The PSF FWHM of 40 arcsec. corresponds to 5 pixels on these 
images. 
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Fig. 5. Images of all the stacked galaxies on the 70 fim COSMOS CLEANed maps with 14 redshift bins. Images are 244 x 244 
sq. arcsec. wide at 70 fim (with 4 arcsec pixel plate). The PSF FWHM of 18 arcsec. corresponds to 4.5 pixel on these images. 



redshift bins: 0<z<0.65 with 317 sources (of which 9 AGNs), 
0.65<z<1.3 with 575 sources (45 AGNs), 1.3<z<2 with 259 
sources (38 AGNs) and z>2 with 198 sources (44 AGNs). 
These bins have been chosen to maximize the number of 
sources present in each bin, while keeping the same Az width. 



2.2. COSMOS data 

The Cosmic Evolution survey (COSMOS) data were acquired 
by MIPS at 24, 70 and 160 fim. The 24 fim observations 
of the COSMOS field is part of two General Observer pro- 
grams (PI D. Sanders): G02 (PID 20070) carried out in January 
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Table 1. Number of 24 //m sources per redshift bins for the COSMOS field. N t0 , is the total number of galaxies used in the stacks, 
and Nacns is the number of sources identified as AGN that were used in the stacks to estimate the AGN contribution. 
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z >2.05 


N M 


2670 


1401 
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1311 


1519 
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2833 


Nagn.s 


83 


76 


55 


225 


232 


200 


288 



2006, G03 (PID 30143) carried out in 2007. We use a to- 
tal net area of 1.93 Sq. degrees. [Le Floc'h et al. (2009)| ex- 
tracted a catalogue at 24 /mi and provided us a sample of 
32840 galaxies with 24 jjm flux densities S 24 > 80 /dy. 
The completeness limit is at the order of 90% at this level. 
(Notice that the survey sensitivities are similar in the COSMOS 
and GOODS fields, at 24, 70 and 160 /mi). The 24 /mi 
galaxies have been completely identified, and redshifts de- 
rived by |Ilbert et al. (2009)1 and |Salvato et al. (2009)| for the 
optically and X-Ray selected sources of the COSMOS field. 



We use the Salvato et al. (2009) photometric redshift catalogue 


of the Cappelluti et al. (2009) 


X-Ray sources catalogue, op- 


tically matched by Brusa et a 


1. (2007)| dBrusa et al., 2010|) to 



identify the AGN in the COSMOS field (Le Floc'h et al., 
2009). Notice that the X-ray flux limits used in the soft 
(0.5-2keV), hard (2-10keV) or ultra-hard (5-10keV) bands 
are 5 x 10~ 16 , 3 x 10~ 15 and 5 x 10~ 15 erg cnr 2 s -1 , 



respectively (|Cappelluti et al., 2007 Cappelluti et al., 2009 



Salv ato et al., 2009| ). We complete this sample with sources 
with a power-law SED ( Alon so-Herrero et al., 20 06) in the red- 
shift range 1.5<z<2.5 using IRAC colors (at lower and higher 
redshifts, the colors can be contaminated by the PAH or stel- 
lar bumps) in the same way as for the GOODS sample. We 
obtained 1668 sources (1115 X-Ray sources, 553 power-law 
sources) detected at 24 /zm and identified as AGNs for 31172 
star forming systems. 

The COSMOS sample being larger than the GOODS one 
we used 14 redshift bins, described in tableQ] The source statis- 
tics in these fields is summarized in Fig.[TJ showing the num- 
ber counts of the GOODS survey (CDFS and HDFN) as well 
as the COSMOS field, corrected for incompleteness. The er- 
rors bars used only include Poisson statistics, and not cosmic 
variance, and are thus likely underestimated. There is no evi- 
dence of a relative major over- or under-density, except maybe 
a slight over-density in the HDFN around 1 mJy, which has a 
negligible contribution to the total background. 



3. Analysis 

3.1. Stacking analysis 

In order to estimate the contribution of mid-infrared galax- 
ies to the 70 /mi and 160 /mi background, we make use of a 



stacking analysif] ( [Dole et al., 2006[ |Bethermin et al., 2010b| l. 
This method consists in stacking the 70 and 160 /mi 
maps at the location of the galaxies detected at 24 /mi. 
The use of this method is justified by two reasons. 1) 
The 24 fim population is a good proxy for the 70 /mi 
and 160 //m populations making up most of the CIB 
near its peak ( |Dole et al., 2006| |Bethermin et al„ 2010] >. 2) 
Only a few sources are individually detected at FIR 
wavelengths and do not resolve much of the background 
dDole et al., 20041 |Frayer et al., 2006a| |Frayer et al., 2006b[ 
|Dole et al., 2006| |Frayer et al., 2009[ >. Notice that stacking 
may suffer from the contamination of galaxy cluster- 
ing, since the stacked image shows in two dimensions 
the 2 point angular correlation function ( [Dole et al., 2006, 
IBavouzet, 2008[ |Bethermin et al., 2010| >. However with the 
Spitzer and Herschel beams, the effects of clustering in the 
stacking are not important (less than 15%) ( IBavou zet, 2008, 
|Fernandez-Conde et al., 2008, Fe rnandez-Conde et al., 2010] >. 

We first stack the 70 /mi and the 160 /mi MIPS data 
(CLEANed maps) as a function of the 24 //m flux, regardless 
the redshift of the sources (Fig.|2]i. This allows us to check the 
consistency of the procedure, since the total brightness mea- 
sured for stacks down to S 24 = 80 //Jy should be equal to the 
sum of the brightnesses obtained by redshift slices, as well as 
identify possible biases. The stacks in the COSMOS and the 
two GOODS fields at 70 and 160 /mi show strong dependencies 
on the fields: while COSMOS and GOODS-N (HDFN) stacks 
are consistent within 20%, stacks in CDFS field is systemat- 
ically lower than COSMOS by a factor of about 1.4-1.8. The 
better statistics in the COSMOS field (surface area and num- 
ber of sources) is limiting the impact of the variance due to 
the large scale structure, and we attribute the systematic lower 
values in GOODS-S to this effect. 

Prior to stacking the 24 /mi catalog in redshift bins on 
the 70 /mi and 160 /mi maps, we use the clean algorithm 
(Hogbom , 1974| > to subtract the few resolved sources present 
in the Far-Infrared maps, in order to remove any bias in the re- 
sulting photometry of the stacked images. The stacking anal- 
ysis presented on Fig. [2] was also done on the cleaned Far- 
Infrared maps. In the COSMOS field, we clean the sources 



1 The IAS stacking library, written in IDL, is publicly avail- 
able at http://www.ias.u-psud.fr/irgalaxies , cf Bavouzet (2008) and 
IBethermin et al. (2010)1 
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Fig. 6. Differential 160 /mi background brightness dvi6o-Bi6oAfe as a function of redshift, in units of nW.rrr 2 .sr _1 , in the COSMOS 
field (left) and GOODS field (right). The lower part of the plot shows a linear zoom between and 0.5 nW.m _2 .sr _1 in order 
to exhibit the trend of cleaned sources and AGNs. Black solid lines: total contribution of infrared galaxies. Solid red lines: 
contribution from resolved sources only. Solid blue line: contribution from AGN only. Dash line Lagach e et al. (2004)| model, 
with S 24 > 80 juJy cut. Dot line:[Le Borgn e et al. (2009)| model, with S 2 4 > 80 /Jy cut. Dot-dash line jValiante et al. ( 2009) 
Dot-dot-dot-dash line: |Befhermin et al. ( 2010c) model, with S24 > 80 /dy cut. Models are discussed in Sect. l4.5l and figure [TUl 



brighter than 80 mJy (resp. 20 mJy) at 160 /mi (resp. 70 /mi), 
level corresponding to 90 to 95% completeness, levels com- 
puted by Monte-Carlo simulations on the data themselves 
(Bet hermin et al., 2010) 1. In both GOODS fields, we removed 
all the detected sources at 160 /mi & 70 //m identified at 24 /mi. 
It corresponds to sources brighter than 19 mJy at 160 /mi (5 
sources in GOODS HDFN & 12 sources in GOODS CDFS) 
and 4.4 mJy at 70 /mi (8 sources in GOODS HDFN & 17 
sources in GOODS CDFS). These brightest detected sources at 
70 and 160 /mi have been individually identified at 24 /mi with- 
out ambiguity, and the redshift of the 24 /mi source is used. The 
flux densities of the removed detected sources are converted 
into brightnesses, and added at the very end of the process to 
account for their CIB contribution (even if it's a small fraction 
at far-IR wavelengths). 

We estimate the AGNs contribution to the CIB as a function 
of redshift using the identifications described in sections 12.11 
andO 

The stacking procedure is performed for each redshift bin 
independently, and the images of the stacks are presented in 
Fig. Hand Fig.g] for the GOODS and COSMOS fields respec- 
tively, together with the measured signal-to-noise ratios. 

We summarize our approach: 

- compute the brightness (by redshift bin) of the detected 
sources that are removed from the maps to create the 
cleaned maps; 

- select galaxies at 24 /mi (all of them, or just AGN, or just 
non-AGN), by redshift bin; 



- stack at the positions of the selected galaxies in the 70 and 
160 fim cleaned maps; 

- perform photometry and bootstrap on those stacks; 

- compute uncertainty budget. 

All the measurements discussed in this section, i.e. num- 
ber of stacked sources, resolved sources, AGN, and resulting 
brightnesses as a function of redshift, are summarized in tables 
iHl and|5] 

3.2. Photometry and uncertainty estimations 

We perform aperture photometry on the stacked images with 
the following parameters at 160 //m: aperture radius of 25", a 
sky annulus to estimate the background between 80 and 110", 
and an aperture correction of 2.29. At 70 /mi, the parameters 
are: 18", 50 and 70", and 1.68. We have secure detections in all 
redshift bins at 160 /mi and 70 /mi, except in the two highest 
redshift bin (z > 1.3) in GOODS. The signal to noise ratio is 
better in COSMOS than in GOODS, due to the larger number 
of sources used (the number of sources used are reported in 
tables [2] to |5}- We thus will take into account in our analysis 
the redshifts bins < z < 1.3 in GOODS and < z < 2.5 in 
COSMOS. 

The error bars come from three quadratic ally summed 
terms. 1- the photometry uncertainty. 2- the Poisson noise com- 
ing from the number of stacked sources. 3- a bootstrap analysis. 

The bootstrap analysis is done by running the stacking pro- 
cess N b times (usually N b = 5000 and N b = 14000 for GOODS 
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Fig. 7. Differential 70 //m background brightness dvjoBjo/dz as a function of redshift, in units of nW.m _2 .sr _1 . in the COSMOS 
field (left) and GOODS field (right). The lower part of the plot shows a linear zoom between and 0.25 nW.m _2 .sr _1 in order 
to exhibit the trend of cleaned sources and AGNs. Black solid lines: total contribution of infrared galaxies. Solid red lines: 
contribution from resolved sources only. Solid blue line: contribution from AGN only. Dash line Lagach e et al. (2004)| model, 
with S 24 > 80 juJy cut. Dot line:[Le Borgn e et al. (2009)| model, with S 2 4 > 80 /Jy cut. Dot-dash line jValiante et al. ( 2009) 
Dot-dot-dot-dash line: |Bethermin et al. ( 2010c) model, with S24 > 80 //Jy cut. Models are discussed in Sect. l4.5l and figure [TUl 



and COSMOS resp.) of a new sample composed of randomly 
selected sources from our original sample, keeping the total 
number of sources constant (Bav ouzet, 2008[ ); this means that 
some stacked positions might be present zero, or multiple times 
in each realization. The bootstrap error bar comes from the 
standard deviation of the distribution of the photometry mea- 
sured on these Nb realizations. Notice that the signal-to-noise 
ratio of the detections in the stacked images (only photomet- 
ric) is higher than the value quoted in this paper, since we add 
the Poisson and bootstrap terms to estimate the final error bar, 
which takes into account the dispersion of the underlying sam- 
ple. The final error bar is thus larger than just the photometric 
noise estimate. The error bars on the AGNs samples were de- 
termined using a smaller number of bootstrap, Nb = 100 and 
N h = 2000 for GOODS and COSMOS respectively. 

The variance due to the large scale structure (also known 
as cosmic variance) and field-to-field variations is a systematic 
component of the noise, that is difficult to estimate at this stage. 
The Poisson noise, used here, gives a strict lower limit of the 
cosmic variance. 



B 



fos = 4.95 + 0.08 nW.m~ z .sr l at 70 /mi (see also the 
summary in Tab.[6]l. 

If we compare with the models from Lagache et al. (2004) 
|Le Borgne et al. (2009)| and |Bethermin et al. (2010c)| (cf sec- 
tion 14.51 and Tab. [6} applying the same selection of using the 
24 /zm sources with S 24 > 80 //Jy, we obtain that we resolve 66 
to 89% of the 160 //m background, and 75 to 98% of the 70 //m 
background. 

To be more specific, using the only post-Herschel model in 
hand (Bether min et al., 2010c| ), our data show that we resolve 
in COSMOS 90% at 160 /mi and 98% at 70 //m of the back- 
ground we are supposed to measure with the selection at 24 /mi 
applied. Our selection introduces an incompleteness in the CIB 
estimate due to the fainter 24 //m sources (S 24 < 80 //Jy), miss- 
ing in our analysis; this loss implies that we resolve 68% of 
the total 160 /zm background and 8 1 % of the total 70 //m back- 
ground in COSMOS (see sect. l4.5l for the details). For compar- 
ison, |BertaetaT(20T0)] identified about 50% of the 100 and 
160 //m backgrounds in individual sources, and account for 50 
to 75% of the background when stacking at the positions of 
24 //m galaxies, as we do. 



3.3. Measurements 

By adding the brightness obtained from the stacking of 24 //m 
sources with S 24 > 80 //Jy and the few detected far-infrared 
sources, we measure B\(,o m _ GOODS = 7.53 + 0.52 nW.m^.sr -1 
at 160 /mi, B 10m _ GOODS = 3.97 + 0.17 nW.m^.sr 1 at 70 



//m and Bi60,„ 



= 7. 



0.19 nW.m ^jr 1 at 160 /mi, 



4. Discussion 

In the following, we present the measurements and the mod- 
els in the form of d(v d ^ versus redshift z, where vB v is the 
CIB brightness in nW.m _2 .sr _1 , A is the wavelength (70 //m 
or 160 //m) and v the corresponding frequency. This represen- 
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tation has the advantage of being independent of the redshift 
binning, thus allowing a direct comparison between datasets 
and models differently sampled in redshift. We discuss data and 
models with the prior selection of S24 > 80 yuJy, and will show 
(sect. 14. 5b that our conclusions for z < 1.5, i.e. where most of 
the FIR background arises, are not modified by this prior selec- 
tion compared to taking fainter galaxies. 

4.1. The 160 fim background: history since z = 2 

The distribution of the 160 fim CIB measured brightness as a 
function of redshift (Fig. [6} shows a plateau between redshifts 
0.3 and 0.9 in COSMOS and GOODS fields, followed by a de- 
crease at higher redshift. The small dip at z — 0.5 in COSMOS 
is not significant, since it disappears when increasing the size 
of the redshift bin (Az = 0.3 instead of 0.15) and is likely due 
to a structure in the COSMOS field. The GOODS field has the 
same trend in redshift. 

The contribution from resolved sources is maximum at z < 
0.3 and strongly decreases at higher redshift, in agreement with 



the identifications of Frayer et al. (2006a)| The AGN contribu- 
tion is rather constant with redshift; the relative contribution of 
AGN thus rises with redshift. Assuming the COSMOS field is 
representative of the whole CIB population, we derive that 33% 
of the 160 /im background is accounted at redshifts < z < 0.5, 
41% for 0.5 < z < 1, 17% for 1 < z < 1.5, and 9% for 
1.5 < z < 2. Our results are consistent with Berta et al. (2010)] 
who analyzed a deep sample in the GOODS-N field at 160 fim 
with PACS/Herschel. Most of the far-infrared sources are re- 
solved by Herschel, and the stacks of 24 fim sources provide 
slightly more depth. Their peak at z — 1 is more pronounced 
than our analysis. 

The |Lagache et al. (2004)| 



Le Borgne et al. (200 9)|Valiante et al. (20 09) and 
Bethermin et al. (2010c) models are overplotted to our mea- 
surements in Fig. [6] using the same selection of S 24 > 80 fi]y 
as applied on the data. Pre-Herschel models predict different 
redshift distribution of dv\^B\a)ldz. The jLagache et al."(20 04) 
models peaks at z ~ 1, and does not fit our data. Our data are in 



better qualitative agreement with the Le Borgne et al. (2009) 
model, (except for z < 0.3), but none models fit the z > 1 
tail. The problem of the discrepancy between our data and the 
Le Borgne et al. (2009) model at z < 0.3 might be twofold: 



our data lacks of statistics at very low redshift due to the 
relative small sky area, and the model might be overpredicting 
low-z galaxies because of the lack of a cold component in the 
galaxies SED used. The Bether min et al. (2010c)| fits well the 
low, intermediate and high redshift ranges, likely because it is 
based on the minimization of recent Spitzer and Herschel data, 
and already takes into account the FIR and submm statistical 
properties of galaxies (see Sect. 14.5b . 
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Fig. 8. The relative contribution of identified AGN to the 70 
and 160 //m CIB as a function of redshift, with the selection at 
24 fim with 5 24 > 80yuJy. 



ments. The peak contribution at 70 fim lies at lower redshift 
than at 160 fim, which is expected as a consequence of the K- 
correction (the effect of the redshifted shape of the galaxies 
spectra). This effect is also seen between 100 fim and 160 fim 
in the PACS/Herschel data by |Berta etal. (2010)| The dip at 
z ~ 0.5 is likely a cosmic variance effect, as it is not seen in the 
GOODS field, and it disappears when we use broader redshift 
bins. This dip does not change our conclusion about the 70 fim 
background emission with redshift. 

The contribution from resolved sources is maximum at 
0.15 < z < 0.3 and strongly decreases at higher redshift, in 



agreement with the identifications of |Frayer et al. (200 6a) The 
AGN contributions is rather constant with redshift. Assuming 
the COSMOS field is representative of the whole CIB popula- 
tion, we derive that 43% of the 70 fim background is accounted 
at redshifts < z < 0.5, 38% for 0.5 < z < 1, 13% for 
1 < z < 1.5, and 5% for 1.5 < z < 2. 



The Lagache et al. (2004) model predicts a peak of the 
CIB at 70 fim at around z - 1, not seen in the data. The 
Le Borgne et al. (2009) predicts a peak at lower redshift (z < 
0.5), with a strong contribution at z ~ 0, not seen either in the 
data; otherwise, the decrease at z > 0.5 has a shape comparable 
to the data, despite a larger high redshift tail. The post-Herschel 
|Bethermin et al. (2010c)| model nicely follows observed evolu- 
tion in COSMOS and GOODS (Fig.©. 



4.2. The 70 nm background: history since z = 2 

The distribution of the 70 fim CIB measured brightness as a 
function of redshift (Fig. [7} shows a maximum contribution 
z < 0.5 in COSMOS, consistent with the GOODS measure- 



4.3. Role of AGN 

The AGN contribution in each redshift bin is shown on Fig. [7] 
and [6] as the lower blue area. The estimate of this contribution 
can be considered as a "best effort" estimate, because of the 
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Fig. 9. The Cosmic Infrared Background history since z ~ 2 at 24, 70, 100 and 160 yum as measured from galaxies selected 
at 24 yum with S 24 > 80/Jy. 24 yum data come from |Le Floc'h et al. (2009)] 70 and 160 pun from this work, all using the 
Spitzer/MIPS data of COSMOS field. 100 and 160 //m data in four redshift bins come from |Berta et al. (2010)| in GOODS-N 
with Herschel/PACS (and using the same 24 yum prior.). 



difficulty of the task of identifying the origin of the far-infrared 
emission (star formation or AGN). Our AGN identification re- 
lies on X-ray detections and IRAC colors (Caputi et al., 2006 
|Salvato et al., 2009), but the far-infrared emission is not neces- 
sarily physically linked to the AGN ( |Le Floc'h et al., 2007) . 

Our analysis shows that the absolute contribution of AGN 
to the CIB at 70 and 160 yum is rather constant with redshift. 
Figure [8] shows the relative contribution to the CIB, and is ob- 
tained by dividing the AGN contribution to the total CIB con- 
tribution. Because of the smaller contribution of higher redshift 
sources to the CIB, the AGN fraction contribution to the CIB is 
increasing with redshift, from about 3 + 10% for < z < 1.5, 
and up to possibly 15-25% for z > 1 .5, but our large error bars 
do not allow any meaningful estimate. We thus can only state 
that the relative AGN contribution is less than about 10% at 



z < 1.5. Our results are in agreement with Dad di et al. (2 007) 
who predict that the contribution of AGNs shouldn't exceed 
more than 7% up to a redshift unity. 



AGN are thought to play a central role in 
terms of physical processes driving galaxy evolu- 
tion and regulating star formation trough feedback 

Ferrarese & Merritt, 2000[ 



(e.g. Magorrian et al., 1998 



|Bower et al., 2006| |Hopkins et al., 2006] |Cattaneo et al., 2009[ 



Hopkins et al., 2010 1. However, our work confirms that, in 
terms of total energy contributions to the CIB, the AGN 
play a minor role. This conclusion is in agreement with the 



identifications of Frayer et al. (2006a)| and the analysis of 
|Valiante et al. (2009)[ where (their figure 19) is shown that less 
than 10% of the sources with S 24 > 80 yt/Jy (i.e. the sources 
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making-up the CIB at 24 to 160 //m, see |Dole etal.,"2 006) 
have a significant AGN contribution. 

4.4. The 24, 70, 100 and 160 backgrounds 

The mid- and far-infrared background buildup at 24, 70, 
100 and 160 /mi as a function of redshift is summa- 
rized in figure [9] and is available onlin^l The 24 /mi data 



the additional constraints of the observed number counts, CIB 
SED 



come from |Le Floc'h et al. (2009) 


and were normalized to 


i 1 — i 

2.86 nW.m .sr ( Bethermin et a 


., 2010|), and the data at 



and 160 //m come from this work. At wavelengths larger than 
60 jt/m, the observed buildup sequence shows an increasing 
contribution from z ~ 0.5 to z ~ 1 with increasing wavelength. 
This behaviour is expected, as consequence of the redshifted 
peak in the galaxy spectral energy distributions, or k-correction 
(e.g. ILagache et al., 2004 Lagache et al., 20051). The 24 /mi 



buildup shows a flatter (or broader) distribution in redshift, with 
a maximum contribution around z ~ 1 ; This mid-infrared dis- 
tribution has a relative z > 1 .5 contribution larger than in the 
far-IR, i.e. the decay slope is smaller at 24 /mi than at 70 /mi 
and larger wavelengths. 

A detailed comparison, however, is still difficult because 
of the cosmic variance. The results from Be rta et al. (2010)| are 
based on GOODS-N, an area about 50 times smaller than used 
here. We furthermore showed that large scale structure is visi- 
ble at z < 0.5 in the COSMOS field. 

The fact that most of the CIB between 100 and 
160 //m is identified as being produced by z < 1 



sources is in line with expectations (e.g. |Lagache et al., 2005 



Beth ermin et al., 20 10c). The 24 /mi background has, how- 
ever, a significant fraction from galaxies at z > 1 (30%, ac- 
cording to |Le Floc'h et al.,~2 009). As expected and observed 
(Marsden et al., 2009) 1, most of the submillimeter background 
is made of sources lying at higher redshifts (z > 1.5). This 
wavelength versus redshift dependence of the background can 
be explained by the sum of SED of galaxies at various red- 
shifts, in particular the peak emission in the far-infrared of the 
reprocessed starlight by dust being redshifted. Thus, the SED 
of the CIB is broader and flatter at far-infrared and submillime- 
ter wavelengths than any individual galaxy SED. 

4.5. The models 

These observations can be confronted to models. We use 
three backwards evolution models developed to fit in- 



frared data: Lagache et al. (2004) Le Borgne et al. (2009) 



IValiante et al. (2009)| and |Bethermin et aL(2 010c) 
among an abundant literature (for the most recent: 
|Franceschini et al„ 2008| |Franceschini et al„ 2009| 

|Pearson & Khan, 2009| |Rowan-Robinson, 2009| >. The main 
feature of the Lagache et al. (2004) model is the use of 



two galaxy populations spectral energy distributions (SED), 
normal/cold and starburst galaxies, parametrized by their 
total infrared luminosity. The local luminosity functions are 
fitted, and the evolution in redshift is applied in order to fit 



and CIB fluctuation level. The |Le Borgne et al. (2009) 
model is based on an automated minimization of the difference 
between the model and selected datasets (local luminosity 
functions, number counts, CIB absolute level) with a given 
SED library ( ]Chary et al„ 200T) . The |Valiante et al. (2009)| 



model introduces scatter in the SED by using Monte-Carlo 
runs within an extended library based on observations from 
the Spitzer archive, and containing starburst and AGN con- 
tributions as a function of IR luminosity. Those three models 
used Spitzer data, and were developed prior to the availability 
of Herschel data. Finally, the Bether min et al. (2010c)| model 
is a fully-parametric approach, automatically fitting the latest 
Spitzer, BLAST & Herschel data with Markov chain Monte 



Carlo (MCMC) ( Dunkley et al., 2 005) method, and using the 
Lagache et al. (2004) SED library of two galaxy populations. 

Figures [6] and Qoverplot these three models in the case of 
a prior selection at 24 /mi (with S 24 > 80 /dy) in order to 
be consistent with the data we are comparing with. Figure [10] 
shows the models, with this prior cut, but also without any cut, 
i.e. all the galaxies. The main differences can be summarized 
as follows: 



- the Lagache et al. (2004) model predicts a contribution of 
infrared galaxies to the 70 and 160 //m backgrounds peak- 
ing at z ~ 1, which is not observed; also the predicted dip 
at z ~ 0.5 is not observed either. 

- the Le Borgne et al. (2009) model overpredicts the galaxy 
contributions at z < 0.3, in disagreement with observations, 
despite the poor statistics of the observations; the origin is 
likely a lack of a cold galaxy component at z ~ 0. The 
general shape of the model at z > 0.3 agrees with the data, 
despite predicted peak and high-redshift (z > 1) tail slightly 
higher than observed. 

- the |Valiante et al. (2009)| model has similar trends as 



http://www.ias.u-psud.fr/irgalaxies/ 



Lagache et al. (2004): almost no contribution at low red- 
shift (in disagreement with data at 70 and 160 /mi), and 
a pronounced peak at z ~ 1, not observed. This models re- 
produces well, however, the z > 1 tail at 160 /mi (but not at 
70 /mi). 

- the |Bethermin et al. (2010c)| model provides a better fit 
to the data, likely because its minimization on re- 
cent Spitzer & Herschel data (Bether min et al., 2 010, 
Oliv er et al., 2010} at far-infrared and submillimeter wave- 
lengths already takes into account the statistical properties 
of galaxies in an empirical way. 

- the selection of S 24 > 80 //Jy galaxies to estimate the back- 
ground buildup with redshift produces an almost flat cut in 
redshift to the brightness (comparison of the two lines for 
each model in fig. \W\ except for |Valiante et al. (2009)| at 
160 //m with larger variations). Thus the peak and the struc- 
ture in redshift observed with the S 24 > 80 /dy cut is not 
much affected by this selection, and our conclusions based 
on this cut can be extended to the whole CIB buildup, at 
least for z < 2.5. However, the S 24 > 80 //Jy cut might 
cause a problem of strong incompleteness at Herschel 
SPIRE wavelengths, made-up by higher-redshift sources 
( [Lagache et al., 2005| |Marsden et al., 20 09); A need to use 
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Fig. 10. Differential 70 (left) and 160 yum (right) background brightness dv v B v /dz as a function of redshift, in units of 
nW.m _2 .sr _1 , as predicted by 4 models: Dash black line: |Lagache et a l. (2004) model; Dot magenta line: |Le Borgne et al. (2009) 
model; Dot-dash grey line: |Valiante et al. (2009)[ Dot-dot-dot-dash blue line: |Bethermin et al. (20 10c)| model. For each model, 
the top line corresponds to all galaxies, and the bottom line to the prior selection (used in this work) of S 24 > 80 /iJy. 



fainter 24 /im flux densities is thus required at submm 
wavelengths. 

All these models predict similar star formation rates and lu- 
minosity function evolutions. Our work put stronger constraints 
on the models, which will have to fine-tune either the galaxies 
SED used or refine the luminosity function evolutions. 



5. Conclusion 



As shown by |Le Floc'h et al. (2009)| and our results, the CIB 
buildup allows to break degeneracies present in the models 
(same predicted number counts and CIB SED, but different 
redshift histories for the luminosity functions for instance). 
Using exquisite Spitzer data on one of the widest and deepest 
fields, we are able to measure that the maximum contribution 
of the 70 fim background (as selected with 24 //m galaxies with 
S 24 > 80 //Jy) occurs at z < 0.5 and at 0.3 < z < 0.9 for the 
160 jum background. 

We measure that 74% of the 160 fj.m background was emit- 
ted at z < 1, and 81% at 70 /mi. We also provided an estimate 
of the AGN contribution to the far-infrared background of less 
than about 10% forz < 1.5. 

The comparison with preliminary Herschel/PACS data on 
GOODS-N from |Bertaet al. (2010)| is in line with our find- 
ings, despite the uncertainties due to large scale structure. The 
consistency of the results confirms that the stacking analysis 
method is a valid approach to estimate the components of the 
far-IR background using prior information on resolved mid-IR 
sources. 



The Lagache et al. (2004) model predictions mainly dis- 
agree with the data, since the peak contribution at z ~ 1 is not 
observed. The Le Borgne et al. (2009) model disagrees with 



the data at low redshift (likely due to the SED used), but suc- 
ceeds in reproducing most of the observed trend, despite an ex- 
cess at z > 1.5. The Bethermin et al. (2010c) model is favored 
by the data. 

Our study, combined with those of |Le Floc'h et al."(2 009) 
and |Berta et al. (2010)| can allow to better constrain the models 
of galaxy evolution, since their predictions can strongly vary 
with redshift, despite good fits of the number counts, luminos- 
ity functions and cosmic infrared background spectral energy 
distribution. 

This study, together with forthcoming works to be done on 
Herschel data, will also help refining the models to compute 
the far-infrared and submillimeter emissivity with redshift, 
needed to compute the optical depth for (hundreds of) TeV 
photons. Since the opacity of the Universe for TeV photons 
depends on the infrared luminosity density along the line of 
sight, the buildup history of the CIB has direct effect on the 
TeV photons propagation. We showed that most (~80%) of 
the far-infrared background is produced at z < 1, in a regime 
where many blazars are observed (e.g. Aharon ian et al., 2 006, 
Albert & Magic Collaboration, 2008). The model predictions 
for TeV obscuration models (e.g. |Mazin & Raue, 2007 
|Franceschini et al., 2008[ 
IKneiske & Dole, 2010[ 



|Stecker & Scully, 2009 



Younger & Hopkins, 2010 



Beth ermin et al., 20 1 0c] > could be disentangled at z < 0.3, 
where the CIB impacts high energy photons the most, by 
comparing their CIB buildup history with our data. 
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Table 2. The CIB brightness by redshift range at 160 //m, in units of nW.m .sr , in the case of the GOODS fields. The 
subscript "stack" refers to the measurement of the signal on the cleaned and stacked image; "sources" to the individually detected 
sources; "AGN" to the sources identified as AGN by Caputi et al. (2006) The total number of sources used in this analysis is 

thus N slack N sources^ 





0< z <0.65 


0.65<z<1.3 


1.3<z<2 


z>2 


N stack 


317 


573 


258 


198 


N 

1 ' sources 


10 


11 


3 


2 


NaGNs 


9 


45 


38 


44 


B 160 mck 
Bl60 source! 
B i60 AGNs 
B 160,,,, 


2.18 ±0.45 
0.57 ±0.18 
0.04 ± 0.07 
2.75 ± 0.46 


2.39 ± 0.54 
0.34 ±0.11 
0.17 ± 0.25 
2.73 ± 0.54 


1.13 ±0.35 
0.06 ± 0.04 
0.32 ± 0.22 
1.19 ±0.35 


0.78 ± 0.27 
0.07 ± 0.05 
0.20 ± 0.19 
0.85 ± 0.27 



Table 3. The CIB brightness by redshift range at 70 fim, in units of nW.ra 2 .sr , in the GOODS fields. Terms are defined in the 
caption of Tab. [2] 





0< z <0.65 


0.65<z<1.3 


1.3<z<2 


z>2 


Nslack 


317 


575 


259 


198 


N 

1 ' sources 


19 


5 


1 





NaGNs 


9 


45 


38 


44 


BjO,,,, 


1.29 ± 0.22 
0.77 ±0.18 
0.05 ± 0.03 
2.06 ± 0.24 


1.25 ±0.28 
0.14 ±0.07 
0.1 ±0.12 
1.4 ±0.28 


0.35 ± 0.15 
0.01 ±0.015 
0.11 ± 0.08 
0.36 ±0.15 


0.15 ±0.12 

0.11 ± 0.12 
0.15 ± 0.12 



Table 4. The CIB brightness by redshift range at 160 fan, in units of nW.m^.sr -1 , in the case of the COSMOS field. The subscript 
"stack" refers to the measurement of the signal on the cleaned and stacked image; "sources" to the individually detected sources; 
"AGN" to the sources identified as AGN by |Salvato et al. (2009)1 The total number of sources used in this analysis is thus 

N stack N sources- 





0<z<0.15 


0.15<z<0.3 


0.3< z <0.45 


0.45< z <0.6 


0.6< z <0.75 


0.75< z <0.9 


0.9<z<1.05 


N stack 


2083 


1559 


2853 


2201 


3225 


3590 


3478 


N 

1 ' sources 


56 


40 


18 


9 


11 


12 


6 


NaGNs 


34 


32 


74 


48 


88 


109 


123 


B 160 mck 


0.43 ± 0.04 


0.56 ± 0.04 


0.87 ± 0.05 


0.63 ± 0.04 


0.76 ± 0.07 


0.75 ± 0.06 


0.71 ±0.06 


Bl60 source! 


0.3 ± 0.001 


0.21 ± 0.002 


0.07 ± 0.001 


0.03 ± 0.001 


0.034 ± 0.001 


0.043 ± 0.001 


0.021 ± 0.001 


B i60 AGNs 


0.02 ± 0.006 


0.01 ± 0.005 


0.03 ± 0.01 


0.02 ± 0.01 


0.02 ± 0.01 


0.01 ±0.01 


0.02 ± 0.01 




0.72 ± 0.04 


0.77 ± 0.04 


0.93 ± 0.05 


0.66 ± 0.04 


0.8 ± 0.07 


0.8 ± 0.06 


0.73 ± 0.06 





1.05<z<1.2 


1.2< z< 1.35 


1.35<z<1.5 


1.5<z<1.65 


1.65<z<1.85 


1.85<z<2.05 


z >2.05 


N slack 


2670 


1401 


2044 


1311 


1519 


2073 


2833 


N 

1 ' sources 


2 





2 


2 





1 


2 


NaGNs 


83 


76 


55 


225 


230 


198 


288 


B 160 stack 


0.52 ± 0.05 


0.33 ± 0.09 


0.36 ± 0.04 


0.26 ± 0.03 


0.22 ± 0.02 


0.27 ± 0.04 


0.47 ± 0.04 




0.005 ± 0.0003 




0.006 ± 0.0003 


0.006 ± 0.0003 




0.003 ± 0.0001 


0.007 ± 0.0003 


B 160,, OWs 


0.02 ± 0.01 


0.02 ± 0.01 


0.01 ± 0.008 


0.04 ± 0.03 


0.04 ± 0.03 


0.05 ± 0.03 


0.04 ± 0.03 


B 160,,,, 


0.53 ± 0.05 


0.33 ± 0.09 


0.37 ± 0.04 


0.27 ± 0.03 


0.22 ± 0.02 


0.28 ± 0.04 


0.47 ± 0.04 
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Table 5. The CIB brightness by redshift range at 70 yum, in units of nW.m 2 .sr 1 , in the COSMOS field. Terms are defined in 
the caption of Tab. [4] 





0<z<0.15 


0.15<z<0.3 


0.3< z <0.45 


0.45< z <0.6 


0.6< z <0.75 


0.75< z <0.9 


0.9<z<1.05 


N stack 


2083 


1559 


2853 


2202 


3225 


3590 


3478 


1 ' sources 


77 


82 


48 


23 


13 


11 


3 


NaGNs 


34 


32 


74 


48 


88 


110 


123 



s 7<w 0.31 + 0.03 0.39 ±0.02 0.6 ± 0.03 0.42 ± 0.03 0.47 ± 0.03 0.41 +0.03 0.39 ± 0.03 

5 7oIL 0.31 +0.001 0.34 ±0.001 0.13 ± 0.0005 0.05 ± 0.0003 0.04 ± 0.0003 0.027 ± 0.0003 0.008 ± 0.0001 

B 10aoNs 0.011 ±0.004 0.008 ±0.003 0.015 ± 0.006 0.009 ± 0.003 0.016 ± 0.005 0.015 ± 0.007 0.015 ± 0.008 

5 7 o * 0.62 ± 0.03 0.73 ± 0.02 0.73 ± 0.03 0.47 ± 0.03 0.51 ±0.03 0.44 ± 0.03 0.4 ± 0.03 







1.05<z<1.2 


1.2<z< 1.35 


1.35< z <1.5 


1.5<z<1.65 


1.65<z<1.85 


1.85<z<2.05 


z >2.05 


N stack 
N 

1 ' sources 
NaGNs 


2670 
2 
83 


1401 
1 

76 


2044 
2 
55 


1311 


225 


1519 


232 


2073 

200 


2833 

288 


Bl0, OIITCe , 

B 70 AGNs 

Bl0 m 


0.25 ± 0.02 
0.004 ± 0.0001 
0.14 ±0.007 
0.26 ± 0.02 


0.18 ±0.01 
0.002 ± 0.0001 
0.011 ±0.005 

0.18 ±0.01 


0.17 ± 0.02 
0.003 ± 0.0001 
0.006 ± 0.004 

0.17 ±0.02 


0.1 ±0.01 

0.28 ± 0.012 
0.1 ± 0.01 


0.09 ± 0.02 

0.026 ± 0.015 
0.09 ± 0.2 


0.09 ± 0.02 

0.022 ± 0.01 
0.09 ± 0.02 


0.16 ±0.02 

0.029 ± 0.018 
0.15 ±0.02 



Table 6. The total CIB brightness at 160 /jm & 70yum for the GOODS & COSMOS fields, in units of nW.m" 2 .^ 1 : lines 1 & 2: our 
estimates; line 3: Betherr mn et al. (2010)| CIB measured value using number counts integration; line 4: Bether min et al. (2010)| 
CIB value with extrapolation of the number counts in power-law; line 5: |Lagache et al. (2004)| CIB model value with the con 
straint: S 24 > 80 /vJy; line 6: Lagache et al. (2004)| CIB model value of the total background; line 7: Le Borgne et al. (2009) 
CIB model value with the constraint: S 24 > 80 /iJy; line 8: Le Borgne et al. (2009) CIB model value of the total background; 
line 9: |Valiante et al. (2009)| CIB model value with the constraint: S 2 4 > 80 /Jy; line 10: |Valiante et al. (2009)| CIB model 
value of the total background; line 11: Betherm in et al. (2010c)| CIB model value with the constraint: S24 > 80 yuJy; line 12: 
Bethermin et al. (2010c) CIB model value of the total background. 



160 /im 70 jjm 



Btot-GOODS 


7.53 ± 0.84 


3.97 ± 0.41 


Btot-COSMOS 


7.88 ±0.19 


4.95 ± 0.08 


$ Bethermin 


9.0 ± 1.1 


5.2 ± 0.4 


" BetherminciBesiimaie 


14.6!" 


6 4 +0 - 7 

u- -0.6 


B modelu^iS 24 > SOfiJy) 


11.91 


5.73 


D 

D model Lagache 


14.87 


6.78 


B 'model UBo rgn S£ '24 > 80^/>') 


9.54 


6.65 


Bmodel LeBorg „ e 


13.57 


8.54 


B mo dei VllKam (S 24 > 80/i/y) 


6.84 


4.27 


D 

D model Va i lante 


16.70 


6.98 


B m odeI Be „, ermi „(S 24 > SOfiJy) 


8.82 


5.02 


f> 'model Bethermin 


11.66 


6.09 



